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There are two independent critical exponents that describe the behavior of systems near their 
critical point. However, at the critical point only the exponent rj, which describes the decay of 
the correlation function, is usually discussed. We emphasize that there is a second independent 
exponent rj that describes the decay of the fourth-order correlation function. The exponent rj is 
related to the exponents determining the behavior of thermodynamic functions near criticality via 
a fluctuation-response equation for the specific heat. We also discuss a scaling law for rj . 

It is well known that the critical behavior of a system near its continuous phase transition^ can be completely 
described by only two independent exponents. 2 This description can be expressed by the scaling behavior of the 
singular part of the free energy density 

fsin e =\t\ 2 - a f±(h/\t\ A ), (t-0±) (1) 

where t = (T — T c )/T c is the reduced temperature, T c is the critical point temperature, and h is the external field that 
is coupled to the order parameter. The scaling functions f±{x) depend on t and h only in the combination x = h/\t\ A . 
[xx note change of notation xx] Because the singular contributions to various thermodynamic quantities are found 
from various derivatives of the free energy in Eq. Q, there are only two independent exponents, say a and A, that 
characterize the critical behavior of the system and thus specify its universality class. For example, the dependence 
on t in Eq. (1) leads to a power-law singularity for the specific heat, C s i ng ~ —d 2 f/dt 2 when h = 0: 

C sing = \tr a c ± (h/\t\ A )=c ± (0)\t\- a , (h = 0). (2) 

At the critical point the thermodynamic derivatives are not useful for characterizing the critical behavior as they 
usually take on infinite values, [xx most of the derivatives are infinite, don't understand this sentence xx] The only 
well known exponent that is defined at the critical point is the anomalous dimension exponent ry, which describes the 
decay of the order parameter correlation function: 

G m (r) = (Sm(r)Sm(O)) = (m(r)m(O)) - (to) 2 - l/r rf - 2+ ", (3) 

where Sm(r) — m(r) — (m) is the order-parameter density fluctuation from the mean at point r. The angular braces 
(• • •} denote an ensemble average, and d is the spatial dimension. Note that G m (r) is a two-point correlation function. 
Although we will use a notation appropriate to magnetic systems, the following discussion also is applicable to other 
systems that have a continuous phase transition, such as fluids at their critical point. 

The search for a second independent exponent leads to the fourth-order correlation function in terms of the order 
parameter density with an exponent denned by£ 

G E (r) = (H(r)H(O)) - (H) 2 ~ l/r d - 2 +"', (4) 

where H(r) is the local Hamiltonian or the interaction energy density. Note that the local interaction energy density 
H{r) (~ m(r) 2 ) is quadratic in the order parameter density m(r). 

The exponents rf and rj are related to the usual exponents via a fluctuation-response equation for the order parameter 
susceptibility 

X = Pjd d rG m (r)~f-^~e^, (5) 
and the energy version of Eq. (|SJ) which relates the thermal susceptibility or the specific heat: 

C/k = (3 2 J d*rG E (r) ~ f ~ £ 2 -"'. (6) 

The correlation length £ characterizes the exponential decay of the correlations, (3 = 1/kT, and k is Boltzmann's 
constant. Equation (JjJJ yields the well-known Fisher scaling identity: 4 



7 = (2 - rj)v, 



(7) 
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where 7 is the linear response susceptibility exponent defined by x ~ \A 7 j an d ^ is the correlation length exponent, 
£ ~ Equation |JBJ yields an energy version of Fisher scalingi^i& 

a={2-rf)u. (8) 

The correlation length exponent v appearing in Eqs. Q and (JSJ) must be the same because of the hyperscaling 
hypothesis, which states that close to criticality the correlation length £ is the only relevant length scale, and thus £ 
is solely responsible for the singular contributions to the thermodynamic quantities. Another important outcome of 
the hyperscaling hypothesis is the Josephson scaling identity, 2 — a — dv, which provides another relation between a 
and v and thus gives a constraint on Eq. (JSJ. If we substitute a = 2 — dv into Eq. @, we obtain another scaling law 
involving rj' : 

i/ = d+2-2. (9) 

We see that the critical exponent rj is determined by the correlation length exponent v and the spatial dimension d. 

To verify the scaling equation © and to test its predictions, we first consider the d = 2 ferromagnetic Ising model 
in the absence of an external field for which exact analytic solutions are available In particular, we are interested 
in the four-point correlation function of the form 

lU4(r) = (SlSzSrSr+l) - {SlS2)(S r S r +l) (10) 

which also can be interpreted as a two-point energy correlation, with S1S2 characterizing the energy of a reference 
bond and S r S r +i the energy of a bond at a distance r from the reference bond, as can be seen in Fig. 1. The Ising 
spins take on the values, S r = ±1. The fourth-order correlation function (r), when summed over r for all distinct 
pairs of bonds is related to the specific heat in agreement with Eq. ©: ^ r wn{r) = de/dfl — 1 + e 2 , where e = (S1S2) 
characterizes the mean energy per bond£ 

From Eq. @ we see that rj' = 2 for the d = 2 Ising model, where we have used the exactly known correlation length 
exponent v = 1. We note that when rj = 2 is substituted into the Eq. ©, we find a logarithmic divergence for the 
specific heat as we should. The asymptotic behavior of w^{r) has been calculated rigorously for the triangular and 
the square Ising lattice and has been found to be identical for the two lattices (apart from an amplitude factor of 1/2) 
and has the usual Ornstein-Zernick formsiS 

e -2r/« 

Wi{r) ^— , (r->oo). (11) 

We see that the predicted value rj = 2 agrees with results from exact analytic calculations. 

For comparison the critical exponents relevant to our discussion for the d = 2 and d = 3 Ising models are summarized 
in Table HJ For d = 3 Eq. predicts rj' = 1.822 (3), where we have used the numerical estimate, v = 0.629 (l). 11 
A direct measurement of the ratio a/v = 2 — rj' can be obtained by a finite size scaling analysis^ of the specific 
heat obtained by Monte Carlo simulations^ The slope of a log-log plot of C versus L provides an estimate for 
the ratio ajv: C ~ t~ a ~ L a l v . A reliable estimate of this ratio for d = 3 is a/v = 0.178 (5). 11 Thus we have 
r/' = 2 — ct/v = 1.822 (5), which is in good agreement with the value predicted by Eq. @. 

In summary, the exponent rj' of the four-point correlation function provides a second independent exponent at the 
critical point. Thus {rf , r/} form an independent exponent pair. The discussion presents a convenient framework for 
introducing the concept of fourth-order correlation functions, which are important in the study of viscous liquids 5 
and glasses^ 
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TABLE I: Comparison of various critical point exponents for the d = 2 and d = 3 Ising models. 



System 


Q 


V 


V 


v' a 


lb 

V 


2D Ising model (exact) 





1 


1/4 


2 


2 


3D Ising model (approx) 


0.112(3) 


0.629(1) 


0.037(1) 


1.822 (3) 


1.822 (5) 



a From Eq. |5J 

'From analyticali™/numericaUi calculations 




FIG. 1: Schematic of four-point spin correlations in the d = 2 Ising model. The line segments represent bonds whose energies 
are determined by the product of nearest-neighbor spins at each end. 



